P-8J 


NASA CASE NO. NPO-18757-1-CU 


PRINT FIG. 


1 


NOTICE 


The invention disclosed in this document resulted from 
research in aeronautical and space activities performed under 
programs of the National Aeronautics and Space Administration. The 
invention is owned by NASA and is, therefore, available for 
licensing in accordance with the NASA Patent Licensing Regulation 
(14 Code of Federal Regulations 1245.2). 

To encourage commercial utilization of NASA-Owned inventions, 
it is NASA policy to grant licenses to commercial concerns. 
Although NASA encourages nonexclusive licensing to promote 
competition and achieve the widest possible utilization, NASA will 
consider the granting of a limited exclusive license, pursuant to 
the NASA Patent Licensing Regulations, when such a license will 
provide the necessary incentive to the licensee to achieve early 
practical application of the invention. 

Address inquiries and all applications for license for this 
invention to NASA Patent Counsel, NASA Management Office-JPL, Mail 
Code 180-801, 4800 Oak Grove Drive, Pasadena, CA 91109-8099. 

Approved NASA forms for application for nonexclusive or 
exclusive license are available from the above address. 

Serial Number: 08/141 . 797 

Filed Date: October 13. 1993 NMO-JPL 


(NAS A-Case-NP0-l 8757- 1-CU) 

OF IMPLEMENTING DIGITAL 
PHASE-LOCKED LOOPS Patent 
Application (NASA. Pasadena 
Office) 89 p 


METH00 


Unci as 


N94-29374 


December 22, 1993 



G3/33 0003660 


Serial No. 08/14 1 , 797 

Filing Dale October 13, 199 3 

/'FOLDOUT FRAME V. • AWARDS ABSTRACT Contract No. NAS7-918 

/ j y Contractor: Caltech/JPL 

Pasadena, CA 91109-8099 

Inventor: Scott A. Stephens JPL Case No. 18757 

Jess Brooks Thomas Jr. NASA Case No. NPO-187 57-1-CU 
Contractor: Jet Propulsion 

Laboratory Dated: August 13 , 1993 


METHOD OF IMPLEMENTING DIGITAL PHASE- LOCKED LOOPS 


The present invention comprises a method of implementation of 
a digital phase-locked loop with loop constants derived from 
roots placed in the s-plane on a root-by-root basis in terms of 
root-specific damping and root specific decay rates. The 
method provides improved flexibility in tailoring high-order 
loop performance and the design is more straightforward and 
understandable for high-order loops. The inventive method 
provides a fully digital formulation, free of analog 
complications. Further, loop bandwidth and damping do not 
change for high-gain loops. 

FIG. 1 is a block diagram of a digital phase-locked loop with 
phase and phase-rate updates; FIG. 2 is a graphical illustration 
of the relation between the conventional second-order loop filter 
parameter r and the controlled-root parameter n z in the 
continuous update limit; FIG. 3 is a graphical illustration of 
relations between third-order loop filter parameters r, k and 
controlled-root parameters x 2 in the continuous update limit 
with the shaded area representing values of r, k where all roots 
are real; FIG. 4 is a graphical illustration of the relationship 
between fourth-order loop filter parameters r, k, a and 
controlled-root parameters , X 2 , n 2 in the continuous update 
limit; FIG. 5 is a graphical illustration of transfer functions 
for loops of order 1 to 4 for various values of n z and X 2 ; 

FIG. 6, comprising FIGs. 6a and 6b, is a graphical illustration 
of normalized loop bandwidth B L T versus root reference-scale 
parameter B L T for a third-order loop without, and with computation 
delay, respectively; FIG. 7, comprising FIGs. 7a and 7b, is a 
graphical illustration of transfer functions for first and second 
order loops with phase and phase-rate feedback and no computation 
delay; FIG. 8, comprising FIGs. 8a and 8b, is a graphical 
illustration of transfer functions for third and fourth order 
loops with phase and phase-rate feedback and no computation 
delay; FIG. 9, comprising FIGs. 9a and 9b, is a graphical 
illustration of transfer functions for first and second order 
loops with phase and phase-rate feedback and computation delay 
of one update interval; FIG. 10, comprising FIGs. 10a and 10b, 
is a graphical illustration of transfer functions for third and 
fourth order loops with phase and phase-rate feedback and 
computation delay of one update interval; FIG. 11, comprising 
FIGs. 11a and lib, is a graphical illustration of root-locus 
plots for a second order loop using discrete-update formulation 
(thick line) and continuous-update approximation (thin line) 
as a function of B L T; FIG. 12, comprising FIGs. 12a and 12b, 
is a graphical illustration of mean time to first cycle slip 
for discrete-update loops with supercritical damping, phase and 
phase-rate feedback, no computation delay and sine phase 
extractor; FIG. 13, comprising FIGs. 13a and 13b, is a 
graphical illustration of quantities describing steady-state 
phase error in discrete-update loops with phase and phase-rate 
_ fe edback and no computation delay; and FIG. 14, comprising 
FIGs. 14a and 14b, is a graphical illustration of quantities 
describing steady-state phase error in discrete-update loops 
with phase and phase-rate feedback and computation delay of 
one update interval. 
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The novelty of the invention resides in providing a method 
for establishing loop filter constants from loop roots in a 
digital approach that avoids analog considerations and thus 
makes DPLL design more straightforward and provides more 
accurately-controlled loop performance. 
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METHOD OF IMPLEMENTING DIGITAL PHASE- LOCKED LOOPS 
ORIGIN OF INVENTION 

The invention described herein was made in the 
performance of work under a NASA contract, and is 
5 subject to the provisions of Public Law 96-517 

(35 USC 202) in which the Contractor has elected not to 
retain title. 

TECHNICAL FIELD 

10 

The present invention comprises a method of 
implementation of a digital phase-locked loop with loop 
constants derived from roots placed in the s-plane on a 
root-by-root basis in terms of root-specific damping 
15 and root specific decay rates. The method provides 
improved flexibility in tailoring high-order loop 
performance and the design is more straightforward and 
understandable for high-order loops. The inventive 
method provides a fully digital formulation, free of 
20 analog complications. Further, loop bandwidth and 
damping do not change for high-gain loops. 



BACKGROUND ART 


It is desired to simplify and improve the design, 
analysis and synthesis of high-order digital phase- 
locked loops; to find a technique for selecting loop 
constants in high-gain loops in a manner that leads 
directly to desired damping and loop noise bandwidth. 

Previous analyses of digital phase-locked loops 
(DPLLs) are based on the traditions of analog loops and 
introduce unnecessary analog considerations such as 
loop-filter time constants and uncontrolled gain 
variations. This reliance on analog tradition makes 
digital-loop analysis unnecessarily cumbersome and 
circuitous and impedes the progress of analysts with 
little analog training. Theory for digital loops can 
be rigorously developed from first principle without 
reference to analog concepts. With an appropriately 
formulated "digital analysis," one discovers that DPLL 
theory and design become more straightforward 
(particularly for third and fourth order loops) and 
that loop performance is more accurately controlled for 
"high gain" loops. 
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Previous analyses have begun with the closed-loop 
equation in the "continuous update" (CU) limit in which 
B^T + 0 , where B L is the loop noise bandwidth and T is 
the loop update interval. For sufficiently small B T 

Ju 

5 (e.g., B l T< 0.02), the CU approximation can provide an 

adequate starting point for loop analysis and design. 
When B L T is increased in this model to larger, "high- 
gain" values, however, loop roots can move away from 
their initial "small B^T" paths in unplanned directions 
10 and the loop can diverge from expected behavior. For a 
loop with discrete update (DU) intervals, a solution to 
the loop equation can be developed in which root 
locations follow predetermined paths as a function of 
B l T. This feature, which is an automatic benefit of a 
15 new controlled-root parameterization of the present 

invention, can provide, for example, supercr itically- 
damped response for all allowed B l T values. 

The following U.S. Patents were found in a search of 
relevant prior art but none is deemed to affect the 
patentability of the inventive method hereof. 


20 
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3,772,600 Natali 
4,426,712 Gorsk i-Popiel 
4,769,816 Hochstadt et al 
4,771,250 Statman et al 
5 4,794,341 Barton et al 

4,82 0,9 93 Cohen et al 
4,821,294 Thomas , Jr. 
4,890,305 Devries 
4,975,930 Shaw 
10 5,016,005 Shaw et al 

5,036,296 Yoshida 
5,073,907 Thomas, Jr. 
5,109,394 Hjerpe et al 
5,122,761 Wischermann 



STATEMENT OF THE INVENTION 


The instant invention comprises an approach that 
reverses the conventional procedure of design for 
obtaining loop filter constants. The instant 
innovation is a method in which loop roots are first 
placed in the s-plane on the basis of new root- 
specific-damping and root-specific-decay-rate 
parameters. Loop constants are then calculated on the 
basis of these roots. 

In the present invention, loop-filter constants are 
specified in terms of new loop parameters. By design, 
each of these parameters has a simple and direct 
physical meaning in terms of a useful loop property: 
loop noise bandwidth, transient decay time or damping. 
For example, a simple choice of parameter values will 
automatically give a loop a particular selected loop 
bandwidth and supercr itically-damped behavior (i.e., 
all roots real, negative, and equal). Thus, the need 
to solve for root location as a function of standard 
loop parameters (e.g., B L , r and k for a third-order 
loop) is eliminated and analysis is simplified. The 
new parameterization is made feasible in a practical 
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sense by the fact that digital loops can usually be 
designed so that they do not suffer significantly from 
the effects of gain variations. That is, variations in 
signal amplitude, due to either gain instability or 
5 signal-power changes, can usually be accounted for by 
using a normalized phase extractor. So comprised, a 
"fully digital" DPLL does not require the analysis or 
precautions necessitated in other DPLL designs by 
potential gain variations. A particularly appealing 
10 benefit of the invention is the elimination of the 
equivocal practice of using the loop parameter r as 
both "damping factor" and an overall gain factor. 

The analysis is extended to fourth-order loops 
15 because of the potential advantages of such loops. In 
some spacecraft applications, loop bandwidth can be set 
to a smaller value for a fourth-order loop than for a 
third. Consequently, lower signal strengths can be 
reliably tracked. Fourth-order DPLLs , unlike their 
20 analog counterparts, are easy to design and implement, 
given the new fully-digital formulation. Accurate 
placement of loop roots results from a simple selection 
of parameters values rather than complicated analog 
circuit design. 



Acquisition in third- and fourth-order loops should 
be carefully crafted. The customary approach for a 
third-order loop has been to acquire first with a 
second-order loop with "wide bandwidth" and then "hand 
over" tracking to a "narrow-bandwidth" third-order 
loop. This approach sacrifices the opportunity of 
directly acquiring weaker signals with the "narrow- 
bandwidth" third-order loop. High-order DPLLs unlike 
similar analog loops, can be easily initialized so that 
they acquire directly, without first acquiring with a 
lower-order loop. Furthermore, if sufficient a priori 
information is supplied, DPLLs will start off tracking 
in-lock, with no transients. The necessary a priori 
information (in the form of a signal phase and its 
derivatives) can be supplied by FFT analysis and/or 
predicts. When a direct acquisition scheme is 
incorporated, fourth-order DPLLs gain even greater 
significance as an option to third-order loops for 
acquiring and tracking weak spacecraft signals. 

To establish a foundation for analysis, a high-level 
description of a DPLL is presented by way of 
background. For loops of first to fourth order, the 
new controlled-root parameterization is used to derive 
a CU-limit solution and to develop a general approach 
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from which numerical, controlled-root solutions to the 
DU loop can be derived. Solutions are given for phase 
and phase-rate feedback, with computation delay for 
closing the loop set to either zero or one update 
5 interval. To tie in with traditional analysis, the new 
loop parameters are related to old loop parameters in 
the CU limit. Loop transfer functions are presented 
for each loop order. The technique for direct 
transient- free acquisition is extended to fourth-order 
loops. Finally, results for two measures of loop 
performance, mean time to first cycle slip and steady- 
state phase error, are described. 


10 
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OBJECTS OF THE INVENTION 

It is therefore a principal object of the present 
invention to provide a method for specifying loop- 
filter constants of digital phase-locked loops in terms 
of loop parameters having direct physical meaning as a 
useful loop property. 

It is another object of the invention to provide a 
method of simplifying the analysis of digital phase- 
locked loops by obviating the need to solve for root 
location as a function of standard loop parameters. 

It is still an additional object of the present 
invention to provide a method for parameterizing 
digital phase-locked loops, wherein loop constants are 
computed on the basis of loop noise bandwidth and 
controlled-root independent parameters related to decay 
times and damping. 
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BRIEF DESCRIPTION OF THE DRAWINGS 


The aforementioned objects and advantages of the 
present invention, as well as additional objects and 
advantages thereof, will be more fully understood 
hereinafter as a result of a detailed description of a 
preferred embodiment when taken in conjunction with the 
following drawings in which: 

FIG. 1 is a block diagram of a digital phase-locked 
loop with phase and phase-rate updates; 

FIG. 2 is a graphical illustration of the relation 
between the conventional second-order loop filter 
parameter r and the controlled-root parameter n 2 in 
the continuous update limit; 

FIG. 3 is a graphical illustration of relations 
between third-order loop filter parameters r, k and 
controlled-root parameters rij X 2 in the continuous 
update limit with the shaded area representing values 
of r, k where all roots are real; 
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FIG. 4 is a graphical illustration of the 
relationship between fourth-order loop filter 
parameters r, k, a and controlled-root parameters n 1 , 

A 2 , n 2 in the continuous update limit; 

FIG, 5 is a graphical illustration of transfer 

functions for loops of order 1 to 4 for various values 

2 

of n and X 2 ; 

FIG. 6, comprising FIGs. 6a and 6b, is a graphical 
illustration of normalized loop bandwidth B L T versus 
root reference-scale parameter B^T for a third-order 
loop without, and with computation delay, respectively; 

FIG. 7, comprising FIGs. 7a and 7b, is a graphical 
illustration of transfer functions for first and second 
order loops with phase and phase-rate feedback and no 
computation delay; 

FIG. 8, comprising FIGs. 8a and 8b, is a graphical 
illustration of transfer functions for third and fourth 
order loops with phase and phase-rate feedback and no 
computation delay; 
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FIG. 9, comprising FIGs. 9a and 9b, is a graphical 
illustration of transfer functions for first and second 
order loops with phase and phase-rate feedback and 
computation delay of one update interval; 

5 

FIG. 10, comprising FIGs. 10a and 10b, is a graphical 
illustration of transfer functions for third and fourth 
order loops with phase and phase-rate feedback and 
computation delay of one update interval; 

10 

FIG. 11, comprising FIGs. 11a and lib, is a graphical 
illustration of root-locus plots for a second order 
loop using discrete-update formulation (thick line) and 
continuous-update approximation (thin line) as a 

15 function of B T; 

Li 

FIG. 12, comprising FIGs. 12a and 12b, is a graphical 
illustration of mean time to first cycle slip for 
discrete- update loops with supercritical damping, phase 
and phase-rate feedback, no computation delay and sine 
phase extractor; 


20 
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FIG. 13, comprising FIGs. 13a and 13b, is a graphical 
illustration of quantities describing steady-state 
phase error in discrete-update loops with phase and 
phase-rate feedback and no computation delay; and 

FIG, 14, comprising FIGs. 14a and 14b, is a graphical 
illustration of quantities describing steady-state 
phase error in discrete-update loops with phase and 
phase-rate feedback and computation delay of one update 


10 


interval 
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1.0 DETAILED DESCRIPTION OF THE INVENTION 


2.0 HIGH-LEVEL DESCRIPTION OF A DPLL 

The block diagram in FIG. 1 shows the basic elements 
of a DPLL. The example DPLL shown in FIG. 1 is based 
5 on "immediate update" of the loop filter and feedback 
of phase as well as phase-rate. Alternate DPLL designs 
might feedback only phase rate and/or have a 
substantial computation delay ("transport lag"). The 
analysis below also treats the case when the 
10 computation delay is equal to one update interval. 

An incoming signal is sampled at a high rate (f 0 ) and 
then counter-rotated sample by sample with model phase 
values generated by a number-controlled oscillator 
15 (NCO) as directed by loop feedback. The resulting 
counter-rotated signal# which should have very low 
frequency, is then accumulated over an update interval 
length T. A phase extractor then processes the 
resulting complex sum to produce a value for residual 
20 phase for the given interval (n th ) . For an ideal phase 
extractor, the n th residual phase is equivalent to the 
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difference of the n th input signal phase <j> n and the n th 
model phase 

fifin = <f>n ^m,n 

5 

with each reference to the center of the sum interval. 
Even though actual residual phase can deviate from this 
linear model due to nonlinearity in the phase extractor 
and/or cycle ambiguities, the theory presented here 
10 will be based on the approximation of Equation (2.1). 
This n th residual phase is immediately passed to the 
loop filter to assist in the prediction of the phase 
rate for the (n + l) c ^ interval. The loop filter 
generates an estimate of phase rate, <J> m>n+1 T, in the 
15 form of phase change per update interval. An estimate 
of the next model phase, the (n + l) th , is projected 
ahead to interval center by adding this estimate phase 
change to the previous n th model phase: 


( 2 . 1 ) 


20 


— 4>m,n 4 * 4>m t n+lT 


( 2 . 2 ) 


The (n + l) th model phase, along with estimated phase 
rate, is used to initialize the phase and rate 
registers of the NCO in a manner that causes the NCO to 
25 produce over the (n + l) th 


interval, phase 
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values characterized by said rate and center- interval 
phase. In this manner, the loop is closed and a new 
value for residual phase is produced to repeat the 
process. 

3.0 LOOP THEORY 

3.1 LOOP FILTER 

As suggested above, residual phase from a given 
interval can either immediately enter the loop filter 
and influence feedback in the very next interval or its 
effect can be delayed to allow more time for 
computation. Two cases are considered here: a) 
"immediate" updates with acceptably small computation 
times and b) updates delayed by one update interval. 
Immediate updates are discussed first. 

The use of residual phase from a given interval to 
update NCO phase in the very next interval can lead to 
a "dead time" during which necessary update 
computations are performed, and during which no 
counter-rotated phasors are accumulated. If the length 
of this "dead time" can be made such a small fraction 


25 
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of an update interval that the SNR loss is acceptably 
small, the immediate update of the loop filter becomes 
feasible. For immediate updates, an N th -order digital 
loop filter, using residual phase values 6<j> up to and 
including the n th sum interval, estimates phase rate 
for the (n + l) th interval according to 

( 3 . 1 ) 

n n t n % j 

4>m,n+lT = K\ 6<f) n + -^2 ^ *** 53 53 ^ ^4 53 53 53 ^ ” 

t’sl »=1 ;=1 »=1 j=l t=l 

where n+i T P^ ase change per update interval T and 
the loop constants K.^ are to be specified below. (An 
N th -order loop has N terms.) 


15 In the case of a relatively long computation time, 

loss of data can often be eliminated by performing 
filter computations in parallel with signal 
accumulation in each update interval, and applying the 
result one interval later than the immediate-update 
20 approach. For a computation delay ("transport lag") of 

one update interval, the loop filter estimates phase 
rate by computing during the n th interval the 
expression 

= Ki stn-1 + * 2 £ Hi + % + K < £EX>*+- 

»=1 y=l »=1 ; = 1 1=1 
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Note that the only change between Equations (3.1) and 
(3.2) is that the n th residual phase is not used in 
computing the (n + l) th phase-rate estimate. 

5 In the following discussion, the loop equation will 

be solved on the basis of immediate updates Equation 

(3.1) , while only results are presented for Equation 

( 3 . 2 ) . 

10 

3.2 CONTINUOUS-UPDATE LIMIT 


3.2.1 LOOP EQUATION 


In many applications, the update interval T is much 
shorter than all other filter time scales, and 
considerable insight may be gained by writing Equation 
(3.1) in the "continuous-update limit", T+0. (One can 
easily show that Equation (3.2) leads to the same 
result in this limit.) To facilitate this, define new 
loop constants k.^ as 


Ki = Ki r 


for* = 1,N 


( 3 . 3 ) 
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10 


so that Equation (3.1) becomes 

* m>n+l ~ -*.**»+ *. f ys*<+ 


(3.4) 


t=s 1 
n 


+^H T T, TS *i +^E r E r E ri ^ + 

‘=1 ; = 1 »'=1 ;=1 k —1 

where estimated phase rate has been rewritten on the 
basis of Equation (2.2) under the assumption that the 
NCO is updated in both phase and phase rate. In the 
limit T + 0 , the first term becomes a derivative and the 
sums become integrals: 

(3.5) 

d f* f* f*' f* f*' t 1 " 

= Ki 6<f> + K 2 J dt 6<f> + k$ J J d?t 6<f> + K4 J J J d?t 6<j> + 


15 


Thus, in this limit, the basic equation governing "NCO 
rate" is the same equation that governs the VCO rate of 
an analog loop, given perfect integrators. 


20 


Solutions of the closed-loop equation can be used on 
the theory of differential equations after substituting 
Equation (2.1) and differentiating N - 1 times: 

d N , <<"-> d»- 3 

dt N 9m+Ki ftN-i 9m + K 3 + ...+ K N (j> m = 

d N-l dN . 2 

+K2 dii 7 ^ +... + KN4 i 


( 3 . 6 ) 
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Solution of this differential equation will give the 
behavior of model phase <f> m in response to input phase <f> . 


5 3.2.2 TRANSFER FUNCTION AND LOOP NOISE BANDWIDTH 

To find the frequency response of the loop, take the 
Fourier transform of both sides of Equation (3.6), and 
utilize the relation 

io d n 6(t) 

m («"')"*{*(»)} ( 3 . 7 ) 

where /{} represents a Fourier transform. This 
produces 

( 3 . 8 ) 

15 ((t 27 ri/)^+/cj(i 27 ri/)^ 1 + K2(i2iti /) n ~ 2 -f- . . . 4- /c yv] $ m (i/) = 

(/Ci(t27Ti/)^’ 1 + /c 2 (t27ri/)^’ 2 + . . . + k n ] $(z/) 

where $ (v) and $ (v ) are the Fourier transforms of 

m 

<P (t) and <j> m (t) . The closed-loop transfer function 

20 H(v) is defined by 4> (v) = H(v)$(v), so we have 

m 

H(u) = + /c 2 (»2 +.. . . + kn 

(t‘2 ni/)* + Ki(i2vv) N ~ l + K 2 (i2m/) N - 2 + . . . + k n 


( 3 . 9 ) 
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Plots of the transfer function for various loop orders 
and constants are shown in FlGs. 7 to 10 . The single- 
sided loop noise bandwidth B L for the closed loop is 
defined by 

5 f°° 

2fl £ = / dt/lflWI 3 (3.101 

J — oo 

This integral is evaluated in Appendix A to find B L as 
a function of kj , k 2 ,... continuous- update limit and 
10 the results are summarized in Table 3-1 for loop fourth 

order . 

3.2.3 PARAMETERIZATION OF LOOP CONSTANTS 

Loop bandwidth B L provides us with one parameter for 
parameterizing the loop constants. Conventional loop 
theory (e.g., 2) specifies additional parameters at 
this point (e.g., r, k,...) to complete the 
parameterization. However, since r, k, etc. cannot be 
given direct physical significance for loops higher 
than second order, they are less useful than more 
carefully chosen parameters. 


20 
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To begin the process of defining new parameters, 
define the higher k's in terms of kj : 

Hi = a { k{ for * = 2,iV (3.11) 

5 No loss of generality is suffered at this step since 

three new parameters replace the three k’s. When this 
equation is substituted for the k's in the expression 
for loop bandwidth, one obtains B L as a function of 
the a's, as shown in Table 3-1. This equation can be 
10 solved for k x in terms of B L and the a's, as shown by 

the first line of equations for each loop order in 
Table 3-2. Based on this expression for k x and 
Equation (3.11), the loop constants can be expressed as 
a function of B L and the a's. 

15 

The variables , will now be parameterized 

in terms of physically significant quantities. To 
begin the search for such parameters, return to 
Equation (3.6). When the roots are unequal (i.e., non- 
20 degenerate) , the general solution to the homogeneous 

equation ( ( t ) = 0) is of the form 

*=1 


(3.12) 
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where each s.^ is a root of the characteristic equation 

+ kis n ~ 1 + K2$ n ~ 2 + ... + *// =0 ( 3 . 13 ) 

and where the a A are amplitudes to be determined by the 

5 initial conditions. We can use characteristics of 

these roots to parameterize k. , k.,..., k , so as to 

12 N 

create a more physically interpretable representation. 


Since the loop constants k. are real, we can 
10 parameterize the roots without loss of generality as 




( 3 . 14 ) 

.} = {-/?i(l±r 7l );-/3 2 (l±r7 2 ); -&(l ± ife ); ...} 


15 


20 


2 

where and 3^ are real. For an odd number of roots, 

. 2 
the unpaired root is not given an n parameter. This 

parameterization takes advantage of the fact that 

complex roots of a poly nominal equation always occur in 

conjugate pair when the polynominal coefficients are 

real. Our goal here is to create parameters which 

dictate the relative behavior of the roots once B L 

has been specified. Thus we can choose one root 

factor, 3i , as the reference decay-rate parameter, 

which is to be determined below as a function of B L , 
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5 


and form new parameters which indicate magnitude 
relative to the first: 

0i = fori >2 (3.15) 

Our root parameterization is now given by: 


(3.16) 

{3i,525 ^3* ^4 5 *®65 • • •} { ""01 (l i *?l)» ~0l ^2(l -t ^7 2 ) » Pi ^(l ^ 3 ) » • • • } 


The overall scale factor 0 , which we choose to be a 
10 positive real number for loop stability, is 

proportional to B L in the CU limit, as can be inferred 
from the specific example in Section 3.2.4. 


2 

The n , which we will refer to as the damping 

15 parameters, contain useful information about each root 

2 

pair. The sign of for a given root pair determines 
the damping behavior for that root pair: 

Vi > 0 two real roots: overdamped 
>7? = 0 two real, equal roots: critically damped 
20 Vi < 0 complex conjugate pair: underdamped/oscillatory 

2 

and the magnitude of is a measure of the separation 
of the root pair: 


£+ — 3- 
+ 3_ 



(3.17) 


25 


The X i 's, which we will refer to as the relative 
magnitude parameters, control the relative magnitudes 
of different root pairs. Furthermore, a negative 
(positive) x i for a given root indicates an 
exponentially growing (decaying) solution to (3.6), 
since is chosen to be a positive real value. Some 
interesting values of A's and n's are: 

all A,- = 1 all real parts of roots equal 

all r/? = 0, A< = 1 all roots real and equal (supercritically damped) 

all 77? = — 1,A,* = 1 “standard” underdamped response 

where "standard" underdamped response for a given root 
corresponds to the response of a 2 nd -order loop with 
r = 2 U = 0.707) . 

To express the k*s as a function of these new root 
parameters, equate each term in Equation (3.13) with 
the like term in the same polynominal factored into its 
roots : 


(s - si)(s - s 2 ) . . . (s ~ -Stf) = 0 


( 3 . 18 ) 
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5 


The k's are then given by 

N 



x 

Kj = + ^ ^ 3{ 3j 

*</ 

N 

«3 = - JZ 

\<j<k 


(3.19) 


(3.21) 


*N = (—1)^(^15233 . . . 3tf) (3.22) 

10 When the root expressions in Equation (3.16) and the k 
expressions in Equation (3.11) are substituted into 
these equations for the k's, one can solve for Bj and 
a's in terms of B L , , and n 2 . Results for the o's 
are given for loops of order 1 to 4 in the second line 
15 of equations for each loop order in Table 3-2. To tie 
in with traditional parameters, the a's are also given 
in terms of r, k, and a in the third line in Table 3-2. 
A graphical mapping between the new parameters and 
traditional parameters is presented in FlGs. 2, 3, and 
20 4 for second-, third- and fourth-order loops, 

respectively. For reference, plots of the transfer 
functions for loops of order 1 to 4 are presented in 
FIG. 5 for various values of x. and n 2 . 
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Since the transient response of a loop is 
characterized by the solution Equation (3.12) to the 
homogeneous equation, knowledge of root locations 

provides a basis for predicting transient behavior and 

2 

5 settling time. Because the n i and values, along 

with the loop bandwidth, completely specify the roots 

by location in the complex plane, loop transient 

response is directly selected at the outset when the 

new loop parameters are chosen. For loops of first to 

10 fourth orders, Table 3-3 presents loop constants for 

two likely implementations: a) standard underdamping, 

where all roots have the same decay time (all X = 1) 

2 

and all -1 and b) supercritical damping, where all 

roots have the same decay time and are critically 

2 

15 damped (all = 0). For comparison. Table 3-3 also 

presents the corresponding conventional parameters. 

3.2.4 DETAILED DERIVATION FOR A FOURTH-ORDER DPLL 

20 

To illustrate the method hereof, a fourth-order loop 
will be treated in detail. As indicated by Equation 
(3.16) , the roots for a fourth-order loop are expressed 
in terms of the controlled-root parameters as 


{*i, *3, *<} = { -0i (1 ± r?i); -0 iA 2 (1 ± rj 7 ) } 


(3.23) 
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Using Equations (3.19) through (3.22) and Equation 
(3.23), we can write: 

/Cl = — (*l + *2 + *3 + * 4 ) 

= (2 + 2\2)0i (3.24) 

K3 — 3l$2 + *1*3 + <1*4 + *3*3 + 32*4 + *3*4 

5 =(4A, + (l-^) + A?(l-,fM (3.25) 

/C3 — “"(*1*2*3 + *1*3*4 + *1*3*4 + *2*3*4) 

= (2A2(1 - Vi) + 2A|(1 - rjl))0f (3.26) 

K 4 = *1*2*3*4 

= Al(l - rj?)(l - nl)e* (3.27) 

Equation (3.24) can be used to express (5j in terms of 

10 k . 

*1 * 

01 = 2T2A7 (3-28) 


This expression and Equation (3.11) can be substituted 
in Equations (3.25) to (3.27) and the three resulting 
equations can be solved for the a's to give 


_ 4Aj + (1 - *?i) + A?(l - r}l) 
(2 + 2Aj) 3 

(3.29) 

_2A,(l-,?)+2A5{l-„i) 

(2 + 2Aj) s 

(3.30) 

n _ Aj(l - i7?)(l - u?) 
(2 + 2Aj)4 

(3.31) 

Analytical evaluation of the loop bandwidth integral in 

Equation (3.10) gives loop bandwidth in 
k 1 s : 

terms of the 

_ K.\K2Kz — K\k\ — K,\k,4 + /Cj/Cs — /C 1 /C 2 /C 4 

— /C 3 /C 4 

4(/Ci/C2«3 — K% — /c}/C 4 ) 


(3.32) 
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5 


Substitution of Equation (3.11) into this expression 
allows a solution for k : in terms of B l and the a's: 


/ci = 4 Bl 


<* 2<*3 — a 3 — 


cczaz — + a\az — cc 2 — < 230:4 


(3.33) 


Note that this equation in combination with Equation 
(3.28) shows that Bj is proportional to B L . 


2 2 

Thus, if n ^ / A 2 , and n 2 are specified, one can 
calculate the a's using Equations (3.20) to (3.31). 
These values of a i can then be substituted in Equation 
(3.33) along with B L to obtain a value for kj. This 
value for k 1 , along with the a's, lead to values for 
k 2 , k 3 , and k 4 through Equation (3.11). As shown in 
Equation (3.3) and listed in Table 3-2, the k's are 
multiplied by an appropriate power of T to obtain the 
K's to be used in the loop filter. In this manner, the 
loop constants for a fourth-order loop (in the limit 

T -> 0) can be calculated in terms of loop bandwidth and 

2 

the controlled-root par ame ters n . and A.. 

11 


20 
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3.3 LOOPS WITH DISCRETE UPDATE INTERVALS 


3.3.1 LOOP EQUATION 


5 


10 


15 


A loop with phase and phase-rate feedback and 
immediate updates will be analyzed in the text while 
Appendix B summarizes differences found in the analysis 
of loops with a computation delay of one update 
interval. Equation (2.1) can be used to rewrite 
Equation (3.1) in terms of the model phase and signal 
phase : 


n n t 

A<£ m , n +i + K\ <f> m ,n + Kz ^ 4>m t { + + . . . — 

1=1 »'=1/=1 
n n j 

K\ 4> n + Ki <j>{ + Ks ^ ^2 4>j + . . . 

»=i »=i y=i 

where the operator a is defined by 


(3.34) 


Ax n — X n ~ i 


( 3 . 35 ) 


and where Equation (2.2) has been used to replace 

estimated phase rate <j> T with A$ , , under the 

c m,n+l m,n+l 

assumption that the NCO is updated with both phase and 
phase rate. To convert Equation (3.34) to a difference 
equation, apply the A operator (N - 1) times: 

A <£m,n+l+ K\A** + -Kj A^ 2 $m,n + . . . + = 

KiA N (f> n + K^A^ *<f>n 4 * . . . 4 * Kf/<f> n 


( 3 . 36 ) 
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where N is the loop order. In analogy with Equation 
(3.6), solution of this difference equation will give 
the behavior of model phase $ in response to signal 
phase <)> . 


3.3.2 TRANSFER FUNCTION 


To find the frequency response of the loop, take the 
10 z-transform of both sides of Equation (3.36) to obtain 

z( 1 - z-‘)"? m (z)+(l - z l ) N -'Ki* m (z) + ... + K„$ m (z) = 

(1 - + ... + K„$(z) 

where $ and $ are the z-transforms of <J> and $ , 

m,n n ' 

15 respectively. To reach this expression, we have used 

the relations 


Z{A" x„} = (l-*-i)" Z{ * n} 

and 


( 3 . 38 ) 


(3.37) 


Z{x n +i} — z Z{x n } 


( 3 . 39 ) 
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where Z{} represents a z-transform. Since the closed- 
loop transfer function is defined by 

$ m (z) = H(z ) $(*) (3.40) 

5 we find that Equation (3.37) yields the expression 



D(z) -(z- 1)" 

D(z ) 


(3.41) 


where 

(3.42) 

D(z) =(z- 1) N + (z - l) w ’ l Ki + z(z - 1 ) n ' 2 K, + z 2 (z - \) n -*K z + . . . + z N ' l K N 


15 


The frequency response of the loop is obtained by 
substi tuting 


z — 


c <2rvT 


(3.43) 


in Equation (3.41) where v is the frequency in H z . 
Plots of the transfer function for various loop 
implementations are presented in FIGs. 7-10. 
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3.3.3 LOOP-NOISE BANDWIDTH 


In analogy with Equation (3.10), the single-sided 
noise bandwidth for the closed loop is defined by 


2B l 


*-^y2r*T) ]2 du 


ft ~ 

■ 4 m 


(3.44) 


10 


which can be rewritten as a contour integral in the 
form 


2BlT = f §(z)H(z- l )z-'dz 


(3.45) 


where the closed path is along the unit circle. This 
integral can be computed on the basis of residues 
15 within the unit circle to obtain B L T as a function of 

the poles of the integrand. For simple poles, the 
integral is given by 

2B t T = £{(*- z i )H(z)H(z-')z - i } (3.46) 

I ' Z—+Zi 

20 where the z are poles of the integrand within the unit 

circle. (For cases with poles of order greater that 
first, the residue evaluation must be appropriately 
mod if ied . ) 
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As seen in Equation (3.41), the poles of H(z) are the 
roots of the polynominal D(z) in the denominator of 
H(z) and therefore satisfy the characteristic equation 

(3.47) 

5 d(2 > = (*- 1 ) w + (^- 1 ) - 1 ) (*_ i) ^-3 jr 3 + 

Let the N roots of this polynominal be z ± so that this 
equation can also be written as 

D{z) = (z - zi)(z - z 2 )(z - z 3 ) . . . (z - z N ) ~ 0 (3.48) 

10 Since the N roots of D(z) must all lie within the unit 
circle if the loop is to be stable (see next 
subsection), all poles of H(z^ ) will lie outside the 
unit circle and will not contribute to the contour 
integral. Furthermore, one can easily show that the 

15 residue for the pole at z = 0 is zero since H(z _1 ) -*■ 0 
as z -»-0. Thus the contour integral can be evaluated 
in a straightforward though algebraically tedious 
fashion on the basis of the N poles of H(z) . The 
resulting expressions for B T as a function of the 

JL 

20 roots are lengthy and uninformative, particularly 
for the higher-order loops, and will not be presented 


here 


3.3.4 LOOP CONSTANTS FROM LOOP ROOTS 


To obtain the relationship between the roots and the 
loop constants, first collect terms according to the 
power of z in Equations (3.47) and (3.48). When the 
coefficients of like powers of z are equated, one 
obtains N equations relating roots and loop constants: 




. - Ks 

( 3 . 49 ) 

YL z<z > 
• </ 

=(")-(">- 


( 3 . 50 ) 

Y2 ****** 

\<j<k 

- (a) 


( 3 . 51 ) 

rp Z\Z<x .. .Ztf 

i Zi 

■p-p*- 

- Ki 

( 3 . 52 ) 

n* 

•l-Ki 


( 3 . 53 ) 


where (£ ) is the binomial coefficient. These N 
equations can be used to solve for each of the N loop 
constants in terms of the N roots, z ± . First solve 
Equation (3.53) for K 1 . The result can then be 
substituted in Equation (3.52) to allow a solution for 
K 2 in terms of the z^ Proceeding sequentially in this 
manner through the rest of the equation, one can obtain 
an expression for each loop constant in term of the 
roots. Thus, if the roots are known, the loop 
constants can be calculated. 
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When the contour integral for B^T is evaluated as a 
function of roots for a given loop order, the result 
can be reduced to a form that contains only the 
functions of z i found on the left-hand side of 
Equations (3.49) to (3.53) . When this form is reached, 
B l T can then be easily expressed as a function of only 
the loop constants. The results are presented in Table 
3-4 for loops of order one to four. 

3.3.5 SOLUTION TO HOMOGENEOUS EQUATION 

Solutions to the homogeneous form of Equation (3.36) 
provide information as to the transient behavior an 
stability of the loop. The general non-degenerate 

solution to the homogeneous equation is given by 

N 

firntti — ^ ^ ( 3 . 54 ) 

* 

where the N amplitudes a^ are to be determined from 

initial conditions and where the N complex numbers z ^ 

are again the roots of the polynominal D(z) in Equation 

(3.47) . (To see that these roots provide solutions to 

the homogeneous equation, substitute <j> = z n in the 

m,n 

left-hand side of Equation (3.36), where the right-side 
set equal to zero, and reduce to the form of Equation 


25 
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(3.47).) Thus, the roots from the homogeneous equation 
are also the poles of the transfer function. In order 
for the loop to be stable, the loop constants must 
be set to values that cause all the roots to fall 
within the unit circle. With a modulus less than 1, a 
root cannot cause the homogeneous solution in Equation 
(3.54) to diverge. 


3.3.6 LOOP PARAMETERIZATION 

Loop parameterization in the case of discrete updates 
interval parallels Section 3.2.3 for the CU limit. 

Loop noise bandwidth B L and the same root-location 
parameters are adopted as independent loop parameters. 

The roots will be parameterized in the form 

{21, * 3 , e -A^,(i±„)r. y , J 55) 

where and n i are the N - 1 independent parameters 
specified in Section 3.2.3, with A^= 1. These 

parameters and "normalized" loop bandwidth B L T will 
comprise the N independent loop parameters needed to 
completely specify the loop. 
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In Equation (3.55), the reference decay-rate 
parameter , which will be represented in the 
normalized, dimensionless form 8^ in the DU case, must 
be determined as a function of these N parameters. In 
5 the CU limit, determination of 8j in terms of B L , 

and , could be carried out explicitly, as shown in 
Section 3.2.3. For DU loops, however, the complexity 
of the equations makes a closed-form solution in the 
general case impractical. Thus, a numerical solution 
0 has been carried out by first selecting a value for B x T 

and the N - 1 independent parameters, and ru, and 
then computing numerical values for the N roots 
through use of Equation (3.55). The resulting z ± 
values can be used to compute the normalized loop 
bandwidth, B L T, as shown in Section 3.3.3, and the loop 
constants, as shown in Section 3.3.4. Repeating the 
process in this fashion on the basis of the same X 
and n i values, one can vary the parameter Bi T 
numerically to obtain B L T and the loop constants as a 
function of B L T. 


20 



In general# B L T increases as Bj T increases from zero 
but can go no higher than a loop-specific maximum 
value. Plots of B L T versus BjT are shown in FIGs. 6a 
and 6b for two supercr itically-damped third-order loops 
with phase and phase-rate update, one with a 
computation delay of zero, the other with a computation 
delay equal to one update interval. In the zero- 
computation delay case, B L T can get no higher than 9.5, 
which corresponds to a 6 1 T value of + °° . In the other 
case, B^T reaches a peak at BjT = -ln(3/4) as b : T 

increases, where B L T « 0.3 is the maximum attainable 
value . 

For a given B L T, therefore, one can find the 
corresponding B X T, if any, and the corresponding loop 
constants on the basis of the analysis outlined above. 
Thus, loop constants can be determined for given 
and n i as a function of B L T. Results are presented in 
Tables 3-5 and 3-6 for loops of order one to four, 
given phase and phase-rate feedback, supercritical 
damping and standard underdamping. Once and the 
loop roots are known, the transfer function can be 
computed on the basis of Equation (3.41). Results are 
plotted in FIGs. 7a and 7b through FIGs. 10a and 10b 
for the two standard loop configurations. 
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Even though a general solution to the DU loop has not 
been obtained, the equations can be expanded as a power 
series in B L T to obtain the higher-order terms relative 
to the CU limit. Results are presented in Table 3-7 
and 3-8 for the loop constants as a power series in B L T 
for the standard loops. 

3.4 EXAMPLE OF STRAYING ROOTS AT LARGE B T VALUES 
IN THE CONTINUOUS-UPDATE APPROXIMATION 

As B^T increases in the CU approximation for loop 
constants (e.g., Table 3-3), the path of the loop roots 
can stray from their original intended course and true 
loop noise bandwidth can exceed the B L parameter value 
used to compute the loop constants. This deviation is 
illustrated in both the sT-plane and z-plane in FlGs. 
11a and lib for a second -order DU loop with phase and 
phase-rate feedback, standard underdamping, and no 
computation delay. The thick straight lines, which are 
based on the exact DU solution presented in preceding 
subsections, show the paths taken by standard- 
underdamped roots as (true) B L T increases from 0 to 
1.4. The thin curved lines, which are produced by the 
CU approximation, follow the DU lines until about 


25 
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B l T =0.1 and then curve toward the real axis. (Note 
that the CU-approximation curves are marked by both 
true B l T values and "parameter" B L T values. True B^T 
is the actual effective noise bandwidth for a DU loop, 

5 while "parameter" B L T is the value used to compute the 

loop constants in Table 3-3.) Where the curves 
separate, the CU approximation starts to diverge from 
standard underdamping and loop transient response 
changes. The inset plot illustrates this divergence in 
0 terms of the corresponding damping parameter n , which 

starts at the intended underdamped value of -1 at 
B l T = 0, increases to a critically-damped value of 0 
at true B l T = 0.8, and then approaches +1 at true 
B l T » 1.2. Thus, loop response at high B l T values will 
5 not match original intended response. In contrast, as 

indicated by the thick lines, the DU exact solution 

2 

maintains standard underdamping (i.e.,n i = -1) for 
allowed values of B^T. 
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4.0 TRAN SI ENT- FREE ACQUISITION WITH DPLLS 


10 


15 


20 


If the signal phase and its time derivatives are 
accurately known at start-up, it is possible to 
initialize the loop sums and loop phase so that the 
loop starts tracking in-lock, with no transients. In 
steady-state tracking, residual phase becomes a 
constant (6<|> = <S<J> ) so that Equation (3.1) or 

X 9S 

Equation (3.2) becomes ( 4 . 1 ) 

&4>n + 1 = Ki S4>»» + K 2 6 <f>i + K 3 d<f>j + K 4 "£ ex >* 

*=1 *=1 y=i t=i y=i k=i 

for up to a fourth-order loop. The quantity n c is the 
computation delay, which is zero or one update interval 
for the present analysis. Note that estimated model 
phase rate, <|> m n + 1 T has been replaced by differenced 
input phase since model phase tracks input phase 
exactly, except for a constant offset, in steady-state 
tracking. Based on this expression, higher-order 
differences become 

n-n 0 n—n c i 

& 2< t>n + 1 = KiS#,, + Kz Y2 6<f>} ( 4 . 2 ) 


»=1 

n-n. 


*=i j'=i 


A 3 $ n + 1 = + K4 ^ 6(f>i 


i=l 


( 4 . 3 ) 

( 4 . 4 ) 
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One can easily show the differences are related to 
input phase values by 

A 0 n +1 = 4>n+l — <kn ( 4 * 5 ) 

A 2 <£ n +i = 0n+ 1 — 2 <j> n + 4 > n - 1 (4*6) 

A^n+l = 0n+ 1 3 <f> n + Z4> n -.\ 4> n —2 (4*"0 

A 4 <^ n +1 = <j>n+ 1 — 4^ n + 6<f>n-l ~ 4 4> n -2 4* 4>n-Z (4*®) 

where each phase value is referenced to interval 
center. We can relate these phase values to 
derivatives of incoming signal phase by means of a 

Taylor expansion whose origin is the center of the n th 
interval (t = t n ) : 

*w =*.+(«- ’ + • • • (4 - 9) 

This expansion shows that phase at the center of the 
nearby intervals is given by 

4>n + 1 = <£(*n + T) = 4>n+ T<t>^ 4* — ^ + 24 ^ + ’ * ' ( 4>1 °) 

4>n = 4>{tr% )=4>n (4.11) 

K-x = *(». - T) = 4*~ + Y+™ ~ T *"' ' + 1 ” ■ • • (4 ' 12) 

K-1 = Wn - 2T) = #» - + 4 - 8 Y^l 3 ’ + is^4 4) ~ • ( 4 - 13 ) 

■An— 3 = #(U - 3T) = - 3 2><»> + 9 y^" ’ " 27 T^ 3) + 1 - • • ■ I 414 ) 
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where T is the separation between intervals ("update 
interval"). By substituting Equations (4.10) through 
(4.14) in Equations (4.5) through (4.8), one can easily 


show the phase differences are related to input phase 
derivatives by 



A*„ +1 = rtf’ + + T*“ 3) + frtf ) + - 

(4.15) 


AV» +l =r^) + g^-) + ... 

(4.16) 

10 

a 3 * „ +I = r 3 4 3) -Y^ 4) + - 

(4.17) 


& 4 <f>n + 1 = T 4 <f>fi 4- . . . 

(4.18) 


By respectively equating Equations (4.15) through 
(4.18) with Equations (4.1) through (4.4), one obtains 
a set of equations whose number is equal to the number 
of unknowns, where the unknowns are <S<t> gs and the loop 
sums. Thus, these unknowns can be expressed in terms 
of the derivatives of input phase. 

To complete initialization of the loop, an estimate 
of starting model phase must be computed. To be exact, 
this estimate must account for steady-state phase 
error. For an arctangent phase extractor, tracking 
error is equal to residual phase (neglecting system- 
noise error and possible cycle ambiguities) . For a 
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sine phase extractor, however, steady-state tracking 
error becomes 


(At - An,*).. = are,»n(2jr t<j>„)/2ir 


( 4 . 19 ) 


where phase is measured in cycles. Thus, model phase 
for the n th interval, after accounting for tracking 
error, is given by 


( 4 . 20 ) 



<f>n ~ 6<f> ts 

<f>n — arc8tn( 2x 


for an arctan extractor 
for a sine extractor 


where phase is measured in cycles. 


The solution for a fourth-order loop will be 

presented in detail. If we assume that time 

derivatives of <t>(t) are negligible above , and 

n 

equate Equations (4.18) and (4.4) , we obtain a value 
for steady-state residual phase for a fourth-order 
loop: 

H.. = y/' (4,21) 


The loop sums are determined in a similar fashion. By 
equating Equations (4.3) and (4.17) , one can solve for 
the single sum: 




( 4 . 22 ) 
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Similarly, the double sum is determined by equating 
Equations (4.2) and (4.16) : 


E E % = > - **.. - ^3 "£' I*i) 


(4.23) 


where the single sum is known from Equation (4.22) . 
Finally, the triple sum is determined by equating 
Equations (4.1) and (4.15) : 


15 


20 


n— n« t j 


T 3 r 3 . T‘ (4.24) 


► • J <i rp~4 rpz m 4 

E EE^‘ = srP’tf > + T^ 21 + T^’ + 


1=1 jssl fc=sl 


n— n 0 


n— n 0 t 


*w .. e «*<-*» E E s ^) 
»=1 »=1 ;'=1 


where the single and double sums are determined by 
means of Equations (4.22) and (4.23) . Thus, if a 
priori values of signal phase and its derivatives are 
available at the start of a track, all of the loop 
variables can be initialized through use of Equations 
(4.20) through (4.24). If the a priori information is 
sufficiently accurate, the loop will start tracking in- 
lock, with no transients. 
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Based on these equations for loop sums and model 

phase, the loop can be initialized as follows. Suppose 

a priori estimates of signal phase and its derivatives 

for a phantom " n th interval" at time t n are available. 

Then Equations (4.22) through (4.24) can provide 

estimates for loop sums at interval "completion", while 

Equation (4.20) provides an estimate of model phase at 

interval center, with steady-state phase error given by 

Equation (4.21). To project ahead to the first 

interval to be processed, the (n + l) th , first estimate 

phase rate according to Equation (4.1) using these 

estimates for loop sums and steady-state phase error. 

Model phase for the (n + l) th interval is then computed 

according to Equation (2.2) on the basis of this 

estimated (n + l) th phase rate and n th a priori signal 

phase. The resulting values for the (n + l) th phase 

and rate are then used to initialize the NCO for the 

(n + l) th interval in the usual fashion. At the 

completion of the (n + l) th interval, sums are updated 

in the standard fashion and substituted in Equation 

(3.1) of Equation (3.2) to provide the (n + 2) th 

estimate of rate. (If n c = 1, the estimated steady- 

state phase error 6<f> is used as residual phase a 

s s 

second time). Equation (2.2) then provides an estimate 
of the (n + 2) th model phase. Loop iteration is normal 
after this point, and steady-state lock should be 


achieved . 



The above analysis has focused on a fourth-order 
loop. Similar analysis can be carried out for loops of 
order lower than fourth, and the results are presented 
in Table 4-1. The preceding derivation assumed a 
signal with polynomial phase so that a steady-state 
phase error would develop. Under less ideal dynamics, 
the above initialization process will not eliminate 
transients but can greatly assist direct acquisition 
with higher-order loops. Similarly, if the derivatives 
of signal phase are known, but phase is not, the loops 
sums can be initialized using Equations (4.22) to 
(4.24) , with initial loop phase arbitrarily set to 
zero. Again, loop acquisition will be greatly 
enhanced. 


5.0 TWO MEASURES OF LOOP PERFORMANCE 

5.1 MEAN TIME TO FIRST CYCLE SLIP 

Simulations have been carried out to determine mean 
time to first cycle slip, ( T lg ) , for loops with phase 
and phase-rate feedback, no computation delay, 
supercritical damping, and a sine phase extractor with 
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perfect amplitude normalization. the tracking-error 
criterion for detecting a cycle slip was 
U - <J> m l > 0.75 cycles. After each slip, the loop was 
reinitialized with perfect a priori so that it would 
5 start off in steady-state lock with no transients. A 

Gaussian random-number generator simulated noise for 
the counter-rotation sums. 

Loops of 2 to 4 order have been simulated on the 
0 basis of the loop constants in Table 3-5. Assumed 

values of B L T ranged between 0.02 and 2.0 and loop SNR 
between 0 and 6 dB. Results are summarized in FlGs. 

12a and 12b, where B L times mean time to first cycle 
slip is plotted versus loop SNR. In terms of cycle 
5 slips, loop performance deteriorates somewhat as loop 

order increases, given a fixed loop SNR. For a given 
loop order and loop SNR, however, cycle-slip 
performance improves as B L T increases, as shown in 
FIGs. 12a and 12b where Bt (T ) is plotted versus B. T 

lSt L 

•0 for 2 nd to 4 th order loops, given a loop SNR of 7 dB. 

j 

With a 3 order loop, for example, FIGs. 13a and 13b 
indicates that B L( T lst ) improves by two orders of 
magnitude when B L T is increased from 0.02 to 0.5. 
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As a test of the simulation software, the cycle-slip 
criteria have been made the same as Viterbi's criteria 
in his exact closed form solution (5) for a first-order 
loop in the CU limit. To within a statistical error of 
about 10%, our results for B L(T lst ) agree with 
Viterbi's theoretical predictions up to SNRl = 4 dB, 
the maximum loop SNR attempted. 

5.2 STEADY-STATE PHASE ERROR 

Loop performance at large values of B L T has also been 

assessed in terms of the steady-state phase error 

(SSPE) . In the CU limit, SSPE is proportional to B T N 

for an N th order loop, as can be derived from Tables 

4-1 and 3-2. For large values of B^T, however, the 

SSPE in a root-controlled loop does not decrease as 

B N as B increases, as illustrated in FIGs. 13a, 13b, 
L L 

14a and 14b. FIGs. 13a, 13b, 14a and 14b plot as a 
function of B L T the dimensionless coefficient required 
to multiply the CU-limit form for SSPE. These plots 
pertain to loops of order 1 to 4, with phase and phase- 
rate feedback, with supercritical damping or standard 
underdamping, and with the indicated computation delay. 
At b l T = 0, the coefficient is equal to the CU-limit 
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coefficient relative to the zero-B L T value is a measure 
of the "excess" SSPE relative to the nominal CU-limit 
values. As FIGs. 13a and 13b indicates, for example, 
the SSPE at B l T =0.5 for a third-order, standard 
underdamped loop, is about four times as large as the 
CU limit would predict. 

FIGs. 13a, 13b, 14a and 14b plot effective loop 

bandwidth as determined from SSPE, where "effective" 

N 

denotes the decrease in bandwidth relative to the B L 
model. In FIGs. 13a and 13b, for example, the 
effective bandwidth is about 0.6 times the actual loop 
bandwidth when B^T =0.5 for a second-order loop. 


SUMMARY 


A first-principles analysis of DPLLs has led to a new 
approach for parameterizing loops. Loop constants are 
computed on the basis of loop noise bandwidth and new 
controlled-root independent parameters that have direct 
physical significance relative to decay times and 
damping. In the continuous-update limit, loop 
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constants of loops of first to fourth order are 
obtained in closed form as a function of these new 
parameters. In a solution for a discrete-update (DU) 
loop, however, complexity of the equations leads to a 
5 numerical approach. The analysis has been applied to 

loops with phase and phase-rate feedback, given either 
zero computation delay or a computation delay equal to 
one update interval. With the new parameterization, 
exact selection of loop bandwidth and damping behavior 
0 can be carried out for high order loops, even when B^T 

is large. 

A method for direct, transient- free acquisition has 
been presented. Given adequate a priori estimates of 
5 phase and its derivatives, steady-state signal clock 

can be obtained directly with third- and fourth-order 
loops without first acquiring with lower-order loops. 

Simulations of loop behavior in terms of mean time to 
first cycle slip have been carried out for loops of 
first to fourth order based on the new 
parameterization, including the larger values of B L T. 
For a given loop bandwidth, loops with larger B L T 
exhibited a considerably better (larger) mean time to 
first cycle slip than those with smaller B L T values. 
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Loop performance has been assessed on the basis of 
steady-state phase error (SSPE) . As B L T increases for 
a given value of B L , SSPE is essentially constant for 
small B l T values (e.g., B L T < 0.02) but increases for 

5 larger values of B L T. Plots provide a measure of this 

"excess" SSPE for large B L T values for loops with phase 


and phase-rate feedback 
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APPENDIX A. LOOP BANDWIDTH FROM LOOP 
CONSTANTS IN THE CONTINUOUS- UPDATE LIMIT 


We wish to solve for loop bandwidth 


10 


15 


2 B L = [ iu\H(v) | 
J — oo 


(-4-1) 


in terms of loop constants k if where the transfer 
function is given by 

(- 4 . 2 ) 

Ki(i2nv) N ~ l + Kj(i2m') N ~ 1 + ... + <» 


B(y) = 77 


(i2xi/) N + Ki(i2iri/) N ~ 1 + Ki(i2xi/) N ~ 2 + ... + Kf/ 


Under the change of variable u = 2irv, we have 

+ K,(iu) N ' 2 + 


( 4 . 3 ) 


H(u ) = Z- 

('■“)" + /c 2 (h,)W-J + — + Kn 

4xB l = du ^(u)! 2 = J du H‘( u)Zr(u) 0 1 - 4 ) 

From Gradshteyn & Ryzhik (3.112) (I.S. Gradshteyn and 

I.M. Ryzhik, Table of Integrals, Series and Products, 
Corrected and Enlarged Edition, Academic Press, New 
York, P. 218 (1980) , we have 

f°° . m n 

loo h N (u)h N (-u) ~ aoA N (A ’ 5) 


' -oo 


20 
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wher e 


15 


M«) = + 6 3U ’("-’> + . . . + 

M“) = w" + «,«"•> + . . . + a N 


5 

and 


ai 

03 

as 

• • • 

0 




ao 

a 2 

a< 


0 



A s — 

0 

• 

ai 

03 

• 


0 




• 

0 

0 

• 

0 


0-N 




b i 

b i 

b 3 


bi sr 




ao 

02 

a< 


0 

10 


m n = 

0 

ai 

a 3 


0 




0 

0 

0 

t • » 

as 


Compar ing 

Equa ti 

ons 

(A. 

3) 

and 

(A. 4 


(AS) 

(A.7) 


(AS) 


(A.9) 


an expression for g N (u): 

y N (u) = [/ci(tu)^” 1 4- K2 (*w) N ’ 2 + ... + «//] [«i(— tu)^" 1 + K2(-iu) N ~ 7 + . . . + «at] 


S 


= E u ^" y) 

y=i 

and an expression for h N (u): 


min(N-jJ-l) 

** + 2 E ( -1 )* K i-k K j 

Jt=i 


+ k 


M“) = (+ i ) N u N + (+i) w - 1 k 1 u*-‘ + (+,-) w - j kj U "-2 + ... 


(A.10) 


(A.ll) 


+ k>n 


20 where the choice of +i rather than -i has been made to 

place the roots of h N (u) in the upper half-plane. 
Finally, by equating like terms in (A. 6) with (A. 10) 



56 


5 


10 


15 


and (A. 7) with (A. 11) , we have the coefficients a^ and 
bj in terms of the loop constants k j : 


y (+i) N ~ J Kj, otherwise. 

min(N-jJ-l) 

*>i = K ) + 2 23 ( _1 ) k K y+» 

fc=l 


(A. 12) 
(A. 13) 


One can now express loop bandwidth as a function of 
loop filter constants by combining (A. 4) and (A. 5) to 
obtain 




i M N 

Aoq A n 


(A.14) 


where the right side is given by Equations (A. 8) and 
(A. 9) with the substitution of Equations (A. 12) , and 
(A. 13) . 


Table 3-1 lists results of Equation (A.14) and 
presents B L as a function of ki for the first four loop 
orders. The table also presents B L as a function of 
the a's after substitution of Equation (3.12). 
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APPENDIX B. ADAPTATION OF ANALYSIS 
TO LOOPS WITH COMPUTATION DELAY OF 
ONE UPDATE INTERVAL 


5 


10 


15 


The preceding DU analysis is based on a loop filter 
with no computation delay. Even though the analogous 
derivation for a loop filter with a computation delay 
of one update interval is in many ways a 
straightforward generalization, there are a few 
differences that are worth mentioning. In analogy with 
Equation (3.34), the loop equation for a loop with a 
computation delay of one update interval and phase and 
phase-rate feedback becomes 


tt—1 n— 1 i 

&<f>m,n+l+Kl + ^2 ^ ^ ^ <£ m ,j + . . . = 

»=l t'=l j= 1 


( 5 . 1 ) 


Ki 4>n—l 


+ k* & + Kz <f>j + . . . 

i=l »= 1 3=1 


20 


Based on this equation, we find in a fashion similar to 
Section 3.3.1 through 3.3.3 that the closed-loop 
transfer function is given by 



D(z) - z(z - 1 )" 
D(z) 


(B.2) 
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where the denominator is defined by 

£>(*) e z(z - 1)" + (z - 1)"' 1 jfl + 2(2 - + 2 J (2 - 1) n - 3 K 3 + . . . + z n -'K n (5.3) 

This polynomial has N + 1 roots in contrast to the 
5 zero-computation-delay case in Equation (3.42) and can 

be represented by 

D(z) = (z - zi)(z - z 2 )(z -Z 3 )...(z- z N+l ) = 0 (£ 4 ) 

Upon expansion of D(z), a constraint equation for the 

st M 

10 extra (N + 1) root comes from the coefficient of z w 

when Equations (B.3) and (B.4) are equated: 

N 

ZN+l ( B .5) 

* 

15 Thus if roots z 1 , z N are designated as 

s t 

selectable roots, the (N + 1) dependent root is 
determined by this equation. As explained in Section 
3.3.6, we are free to place the N selectable roots, as 
before, according to loop noise bandwidth, damping 
20 parameter and relative-magnitude parameter. The DU 

analysis of a computation-delay loop is otherwise 
parallel to the analysis for zero-computation-delay 
loop. The linear equations relating the loop constants 
to the independent roots are similar to 
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Equations (3.49) through (3.53) and will not be shown 
here. Results for this implementation have been 
presented in the text. 

5 Having thus disclosed a preferred mode of carrying 

out the invention, what is claimed is: 



Table 3-1. Loop Bandwidth from Loop Constants in the Continuous-update Limit 


l at order 
Bl = ^ 


2 nd order 

Bl = f (1 + «») 


3 rd order 

d _ /c?/c 2 - /ci/c 3 + /c^ v _ /ci 02 - 03 + o£ 
L 4 (/ci/c 2 - /c 3 ) ' 4 « 2 - «3 


4 tA order 

_ *1*2*3 — /Ci/Ca — /C?/C4 + /C3/C3 — /Cf/Cg/C^ — /C3/C4 V 
4 (/Ci/C 2 /C 3 — /C3 — /CJ/C4) ' 


/ci a 2 Q3 - Qa ~ <*4 + tt%a 3 ~ 02Q4 ~ 03Q4 

4 a 2 a3 — a| — at 



Table 3-2. Loop Filter Constants in the Continuous-update Limit 


1 4< order . 


K\ = 4B l T 


2 nd order 


Controlled-root: 

Standard: 


Z rd order 


K 2 = a 2 K\ 



K x = ABjJl — - ° 2 


OL 2 ~ OE3 + <* 2 

Controlled-root: 

Standard: 


4 th order 


K 2 = a 2 K\ 

a 2 - 2A 2 ± (l-^i) 

(2 + A 2 ) 

a 2 = F 


K 3 = a 3 K\ 

a 3 = -^ 

r 


Kx = 4BlT 


oc 2 a 3 - a, - a 4 

x — “ x 

a 2 a3 - a 3 - a 4 - a 2 a 4 - a 3 a 4 + 0^03 
Controlled-root: 


if 2 = a 2 K\ 


K 3 = aaiCf 


if 4 = a^J* 


a . = 4A 2 + (l-i7?)+Aj(l-r?|) _ 2A 2 (l-f7 ?) + 2A|(1-t?|) _ A^(l-r?f)(l-^) 

2 (2 + 2A 2 ) 2 3 “ (2 + 2A 2 ) 4 a4_ (2 + 2A a ) 4 


a 2 = i 


<*3 = \ 

r 


a 4 = -% 

r 


Standard: 



Table 3-3. Loop-filter Constants for Typical Implementations 
in the Continuous-update Approximation 

Supercritically damped: t;? = 0, A, = 1, for all roots 

Loop constants Conventional par am eters 



K x 

K, 

K, 

K, 

r 

k 

a 




(« 3 Kf) 


(1/0:2) 


(a 4 r : 

1 st Order 

4 B l T 







2 nd Order 

fB L T 

-K 2 

4 a x 



4 



3 rd Order 

jjB l T 


h K -- 


3 

1 

3 


4 tA Order 

256 p r r> 

ecloo 

h K 1 

255^1 

8 

5 

4 

5 

2 

27 



Table 3-3. continped Loop-filter Constants for Typical Implementations 
in the Continuous-update Approximation 


Standard underdamped: 

*1i - -It 

A i = 1, for all roots 





Loop constants 
K! K 2 

?) 

K s 

(asKl) 

K< 

(«<«?) 

Conventional parameters 
r k 

(l/o 2 ) (o: 3 r 2 ) 

(Ur 3 ) 

1** Order 

4 B l T 







2 nd Order 

\BlT 

\K\ 



2 



3 rd Order 

§b l t 


$r K f 


9 

i 

3 

5 


4 th Order 

%B l T 




2 

1 

5 

i 



Table 3-4. Loop Bandwidth from Loop Constants for Discrete-update Loops with 
Phase and Phase-Rate Feedback and No Computation Delay 


1** order 


2 nd order 

j3 rp _ 2 If? + 2K 2 -f 
L ~ 2Ki(4 - 2Ki - K 2 ) 


3 rd order 


B l T = 


4K±K2 - jKiKs ±4Kl ± 2KiKl ± 4K?K 3 + 4K 7 K 3 ± ZKiKrKa ±Kl± K x Kl 
2 (K x K 2 - K 3 + K x K 3 )(8- 4K\ - 2Ki - K 3 ) 


4 th order 


/ 8KiKtK a - SKiKi - 8 K\K a + 8 K\K Z - 8K X K 2 K A - 8K Z K A + 4K X K^K 3 + 8KlKl + 8K 2 Kl + 4K\K 2 K a + 8 K\K a - 
2 0K t K z K A - 8K\ + 8K x K 2 Kl + 2 K| + 14K\K z K a + UK 2 K Z K A + 4K X K$K A - lOKtfl + 11K X K 2 K 3 K 4 + 

V 2KiKl +J>KIKa -f- 7K?K[ + 6K 2 K* ± SK X K^K 4 + 5 K x K z Kl + 4 KzKj + 4K x K 3 Kl ± Kl + K X K\ 

~~ 2(K X K 2 K 3 - K'i - K'lK A + K X K£ + K X K 2 K 4 - 2K 3 K 4 + 2K X K 3 K 4 - K* + Ki^)(16 - 8^ - 4K 2 - 2 K z - K 4 ) 



Table 3-5. Loop Filter Constants for a Discrete-update Loop with Phase 
and Phase-rate Feedback and Supercritically-damped Response 
No computation delay 



1st order 

2nd order 

3rd order 


4th order 



b l t 

K x 

K x 

K 2 


*2 


K x 

k 2 

*3 

K< 

0.001 

0.003992 

0.003193 

2.553e-06 

0.002903 

2.812e-06 

9.084€-10 

0.002747 

2.833e-06 

1.299e-09 

2.234e~13 

0.005 

0.01980 

0.01582 

6.309e-05 

0.01438 

6.941e-05 

1.118e-07 

0.01361 

6.988e-05 

1.597e-07 

1.369e~10 

0.01 

0.03922 

0.03130 

0.0002488 

0.02845 

0.0002733 8.778e-07 

0.02692 

0.000275 

1.251e-06 

2.138e-09 

0.02 

0.07692 

0.06125 

0.0009677 

0.05567 

0.001060 

6.765e-06 

0.05269 

0.001065 

9.617e-06 

3.265e-08 

0.03 

0.1132 

0.08993 

0.002118 

0.08174 

0.002312 

2,22Ge-05 

0.07738 

0.002321 

0.0000312 1.578e-07 

0.05 

0.1818 

0.1438 

0.005576 

0.1307 

0.00605 

9.485e-05 

0.1237 

0.006059 

0.0001337 1.113e-06 

0.075 

0.2609 

0.2051 

0.01176 

0.1864 

0.01267 

0.0002936 

0.1766 

0.01265 

0.0004113 5.055c-06 

0.1 

0.3333 

0.2607 

0.01965 

0.2369 

0.02101 

0.0006405 

0.2245 

0.02094 

0.0008915 1.439«-05 

0.15 

0.4615 

0.3572 

0.03931 

0.3248 

0.04147 

0.001847 

0.3080 

0.04113 

0.002540 

5.978e-05 

0.2 

0.5714 

0.4379 

0.06264 

0.3983 

0.06523 

0.00378 

0.3780 

0.06443 

0.005139 

0.0001670 

0.25 

0.6667 

0.5061 

0.08835 

0.4606 

0.09089 

0.006432 

0.4375 

0.08942 

0.008652 

0.0003222 

0.3 

0.7500 

0.5643 

0.1155 

0.5139 

0.1175 

0.00976 

0.4885 

0.1152 

0.01300 

0.0005671 

0.35 

0.8235 

0.6142 

0.1436 

0.5598 

0.1444 

0.01371 

0.5327 

0.1411 

0.01808 

0.0008996 

0.4 

0.8889 

0.6575 

0.1720 

0.5998 

0.1712 

0.01821 

0.5712 

0.1668 

0.02380 

0.001325 

0.45 

0.9474 

0.6952 

0.2006 

0.6348 

0.1977 

0.02321 

0.6050 

0.1920 

0.03006 

0.001844 

0.5 

1.0 

0.7282 

0.2291 

0.6657 

0.2235 

0.02864 

0.6349 

0.2166 

0.03679 

0.002459 

0.6 


0.7828 

0.2851 

0.7173 

0.2732 

0.04059 

0.6852 

0.2634 

0.05133 

0.003967 

0.7 


0.8257 

0.3394 

0.7585 

0.3196 

0.05371 

0.7258 

0.3069 

0.06692 

0.005832 

0.8 


0.8599 

0.3915 

0.7920 

0.3629 

0.06769 

0.7590 

0.3471 

0.08318 

0.008026 

0.9 


0.8874 

0.4414 

0.8196 

0.4030 

0.0823 

0.7865 

0.3842 

0.09983 

0.01052 

1. 


0.9096 

0.4890 

0.8426 

0.4402 

0.09735 

0.8097 

0.4184 

0.1167 

0.01328 

1.2 


0.9425 

0.5779 

0.8782 

0.5065 

0.1283 

0.8462 

0.4790 

0.1503 

0.01951 

1.4 


0.9645 

0.6588 

0.9042 

0.5636 

0.1597 

0.8734 

0.5308 

0.1832 

0.02650 

1.6 


0.9793 

0.7328 

0.9237 

0.6130 

0.191 

0.8942 

0.5754 

0.2151 

0.03410 

1.8 


0.9889 

0.8006 

0.9386 

0.6559 

0.222 

0.9106 

0.6141 

0.2457 

0.04216 



Table 3-5. CONTINUED Loop 

Filter Constants for a Discrete-update Loop with Phase 


and Phase-rate Feedback and Supercritically-damped 

Response 

No computation delay 









2. 


0.9950 

0.8631 

0.9502 

0.6935 

0.2525 

0.9236 

0.6479 

0.2749 

0.05058 

2.5 


1.0 

1.0 

0.9697 

0.7689 

0.3258 

0.9466 

0.7158 

0.3419 

0.07270 

3. 




0.9811 

0.8248 

0.3947 

0.9612 

0.7667 

0.4011 

0.09566 

3.5 




0.9881 

0.8672 

0.4591 

0.9710 

0.8058 

0.4533 

0.1189 

4. 




0.9925 

0.8997 

0.5194 

0.9778 

0.8366 

0.4996 

0.1422 

4.5 




0.9953 

0.9248 

0.576 

0.9827 

0.8613 

0.5409 

0.1652 

5. 




0.9971 

0.9444 

0.6291 

0.9864 

0.8814 

0.5779 

0.1879 

Computation delay = 

1 update 

interval 








1st order 

2nd order 

3rd order 


4th order 



B l T 

K ! 

K ! 

K* 

K ! 

K? 

k 3 

K i 

K 9 

K s 

K 4 

0.001 

0.003976 

0.003181 2.538e-06 

0.002892 

2.796e-06 

9.015e-10 

0.002737 

2.816e-06 

1.288e-09 

2.211&.13 

0.005 

0.01942 

0.01554 

6.135e-05 

0.01414 

6.749e-05 

1.077e-07 

0.01339 

6.795e-05 

1.537e-07 

1.306e»10 

0.01 

0.03779 

0.03023 

0.0002357 

0.02752 

0.0002588 8.165e-07 

0.02606 

0.0002604 1.163e-06 

1.952e-09 

0.02 

0.07177 

0.05734 

0.0008737 

0.05224 

0.0009552 5.897e-06 

0.04951 

0.0009599 8.364e-06 

2.745e-08 

0.03 

0.1027 

0.08185 

0.001832 

0.07461 

0.001993 

1.809e-05 

0.07076 

0.002000 

2.554e-05 

1.231e-07 

0.05 

0.1571 

0.1245 

0.004476 

0.1136 

0.00482 

7.032e-05 

0.1078 

0.004819 

9.836e-05 

7.611e-07 

0.075 

0.2148 

0.1685 

0.008742 

0.1538 

0.009282 

0.0001955 

0.1462 

0.009237 

0.0002700 

3.007e-06 

0.1 


0.2046 

0.01371 

0.1869 

0.01433 

0.0003895 

0.1778 

0.01420 

0.0005309 

7.609e-06 

0.15 


0.2594 

0.02487 

0.2377 

0.02515 

0.0009724 

0.2265 

0.02467 

0.001289 

2.610e-05 

0.2 




0.2740 

0.03595 

0.001784 

0.2618 

0.03495 

0.002295 

5.906«-05 

0.25 




0.3000 

0.04617 

0.002793 

0.2878 

0.04452 

0.003473 

0.0001076 

0.3 







0.3071 

0.05316 

0.004765 

0.0001722 

0.35 







0.3211 

0.06074 

0.006129 

0.0002543 


Table 3-6. Loop Filter Constants for a Discrete-update Loop with Phase 
and Phase-rate Feedback and Standard-underdamped Response 


No computation delay 



1st order 

2nd order 

3rd order 


4th order 



b l t 

K ! 

K i 

k 7 

K x 

Ki 

k 3 

Kt 

k 2 

K 3 

K 4 

0.001 

0.003992 

0.002661 

3.545e-06 

0.002603 

3.016e-06 

1.311e-09 

0.002367 

2.804e-06 

1.661e-09 

4.923e-13 

0.005 

0.01980 

0.01319 

8.752e-05 

0.01290 

7.437e-05 

1.611e-07 

0.01173 

6.907e-05 

2.037e-07 

3.008e-10 

0.01 

0.03922 

0.02609 

0.0003448 

0.02552 

0.0002926 1.264e-06 

0.02321 

0.0002716 1.594e-06 

4.693e-09 

0.02 

0.07692 

0.05106 

0.001338 

0.04996 

0.001133 

9.720e-06 

0.04545 

0.001051 

1.222e-05 

7.147e-08 

0.03 

0.1132 

0.07499 

0.002922 

0.07338 

0.002469 

3.15645-05 

0.06679 

0.002289 

3.955e-05 

S.447e-07 

0.05 

0.1818 

0.1199 

0.007658 

0.1174 

0.006445 

0.0001355 

0.1070 

0.005966 

0.0001687 2.420e-06 

0.075 

0.2609 

0.1713 

0.01607 

0.1677 

0.01345 

0.0004177 

0.1530 

0.01243 

0.0005159 

1.093e-05 

0.1 

0.3333 

0.2179 

0.02670 

0.2133 

0.02226 

0.0009073 

0.1949 

0.02054 

0.001112 

3.096e-05 

0.15 

0.4615 

0.2991 

0.05288 

0.2929 

0.04370 

0.002595 

0.2683 

0.04022 

0.003136 

0.0001273 

0.2 

0.5714 

0.3675 

0.08345 

0.3599 

0.06842 

0.005269 

0.3305 

0.06282 

0.006285 

0.0003313 

0.25 

0.6667 

0.4258 

0.1166 

0.4171 

0.09491 

0.008901 

0.3838 

0.08697 

0.01049 

0.0006738 

0.3 

0.7500 

0.4760 

0.1511 

0.4662 

0.1222 

0.01341 

0.4299 

0.1118 

0.01562 

0.001176 

0.35 

0.8235 

0.5196 

0.1862 

0.5089 

0.1496 

0.01871 

0.4701 

0.1367 

0.02156 

0.001851 

0.4 

0.8889 

0.5577 

0.2213 

0.5463 

0.1768 

0.02470 

0.5056 

0.1613 

0.02817 

0.002705 

0.45 

0.9474 

0.5914 

0.2561 

0.5793 

0.2034 

0.03129 

0.5370 

0.1853 

0.03534 

0.003738 

0.5 


0.6214 

0.2902 

0.6085 

0.2294 

0.03838 

0.5650 

0.2087 

0.04296 

0.00495 

0.6 


0.6721 

0.3560 

0.6580 

0.2788 

0.05381 

0.6128 

0.2532 

0.05920 

0.007884 

0.7 


0.7134 

0.4181 

0.6983 

0.3247 

0.07047 

0.6521 

0.2944 

0.07632 

0.01145 

0.8 


0.7475 

0.4763 

0.7315 

0.3671 

0.08796 

0.6848 

0.3325 

0.0939 

0.01558 

0.9 


0.7762 

0.5306 

0.7593 

0.4061 

0.1060 

0.7124 

0.3675 

0.1116 

0.02021 

1 . 


0.8007 

0.5813 

0.7829 

0.4421 

0.1243 

0.7360 

0.3997 

0.1293 

0.02527 

1.2 


0.8399 

0.6727 

0.8205 

0.5059 

0.1612 

0.7743 

0.4569 

0.1639 

0.03644 

1.4 


0.8699 

0.7524 

0.8490 

0.5603 

0.1979 

0.8038 

0.5058 

0.1971 

0.04869 

1.6 


0.8936 

0.8223 

0.8713 

0.6072 

0.2337 

0.8272 

0.5479 

0.2286 

0.06171 

1.8 


0.9127 

0.8839 

0.8891 

0.6478 

0.2684 

0.8461 

0.5846 

0.2582 

0.07525 


Table 3-6. continued Loop Filter Constants for a Discrete-update Loop with Phase 
and Phase-rate Feedback and Standard-underdamped Response 


No computation delay 


2. 


0.9286 

0.9385 

0.9035 

0.6832 

0.3019 

0.8618 

0.6167 

0.2860 

0.08913 

2.5 


0.9587 

1.051 

0.9297 

0.7544 

0.3802 

0.8910 

0.6817 

0.3482 

0.1245 

3. 


0.9827 

1.134 

0.9470 

0.8076 

0.4510 

0.9111 

0.7310 

0.4013 

0.1599 

3.5 




0.9592 

0.8485 

0.5151 

0.9256 

0.7695 

0.4470 

0.1945 

4. 




0.9680 

0.8805 

0.5735 

0.9366 

0.8004 

0.4867 

0.2281 

4.5 




0.9747 

0.9061 

0.6269 

0.9452 

0.8256 

0.5213 

0.2605 

5. 




0.9798 

0.9267 

0.6760 

0.9520 

0.8465 

0.5519 

0.2917 

Computation delay = 

1 update 

interval 








1st order 

2nd order 

3rd order 


4th order 



b l t 

K ! 


k 2 

Ki 

K a 


K i 

K 2 

Ki 

k a 

0.001 

0.003976 

0.00265 

S.516e-06 

0.002594 

2.997e-06 

1.299e-09 

0.002358 

2.784«-06 

1.645e*09 

4.859e-13 

0.005 

0.01942 

0.01299 

8.487e-05 

0.01270 

7.229e-05 

l.S49«-07 

0.01155 

6.717e-05 

1.959e-07 

2.862e-10 

0.01 

0.03779 

0.02533 

0.0003248 

0.02474 

0.0002769 1.170e-06 

0.02252 

0.0002573 

1.478e-06 

4.259«-09 

0.02 

0.07177 

0.04827 

0.001194 

0.04709 

0.001019 

8.391e-06 

0.04293 

0.0009475 

1.058«-05 

5.935e-08 

0.03 

0.1027 

0.06919 

0.002479 

0.06739 

0.002121 

2.554«-05 

0.06154 

0.001971 

3.210d-05 

2.636«-07 

0.05 

0.1571 

0.1060 

0.005937 

0.1030 

0.005094 

9.769e-05 

0.09425 

0.004732 

0.0001220 

1.594e-06 

0.075 

0.2148 

0.1448 

0.01130 

0.1401 

0.009726 

0.0002657 

0.1285 

0.009020 

0.0003292 

6.120«-06 

0.1 


0.1775 

0.01725 

0.1709 

0.01489 

0.000518 

0.1571 

0.01379 

0.0006363 

1.503«*05 

0.15 


0.2300 

0.02963 

0.2193 

0.02572 

0.001236 

0.2021 

0.02372 

0.001494 

4.860d-05 

0.2 

• 

0.2713 

0.04155 

0.2554 

0.03630 

0.002167 

0.2360 

0.03335 

0.002579 

0.0001037 

0.25 


0.3071 

0.05222 

0.2833 

0.04616 

0.003232 

0.2622 

0.04227 

0.003795 

0.0001781 

0.3 




0.3055 

0.05518 

0.004373 

0.2832 

0.05037 

0.005075 

0.0002685 

0.35 




0.3242 

0.06338 

0.005530 

0.3002 

0.05766 

0.006374 

0.0003710 

0.4 







0.3143 

0.06420 

0.007664 

0.0004821 

0.45 







0.3263 

0.07006 

0.008927 

0.0005986 



Table 3-7. Loop Filter Constants as Power Series in B^T in the Case of Phase and 

Phase-rate Feedback and No Computation Delay 


Supercritically damped: rft = 0, A< = 1, for all roots 

Kt 



n _ Ks 

as = K f 



1*‘ Order 4 B L T - 8(B L T) 2 + 16(£ L T) 3 

2 nd Order $B L T- *g(B L T) 2 + ^(B^T) 3 J + § B L T - $- h (B L T) 2 

3' d Order + § + gf**7- fgg&UWV 

4“ Order »**7 " 8W(***V + WSWUMT ! + 8**7 ~ SSStW 


£ + &B*T+£!gj|(B L T) 9 <* 

VO 

A + ^r+iSyS'- 7 ') 1 5fe + AB t T + 5 S(i> t r)' 


Standard underdamped: 17? = -1, A< = 1, for all roots 


K , 


K, 

02 = 4 


aa = 4 


a4 = 4 


1“ Order 4B L T - 8[B L T) 2 + 16{B L T) 9 
2 nd Order \B L T - $?(. B L T ) 2 + ^(BiT) 3 
3 rd Order I %B L T- %jg{B L T)* + ^S(B^) 3 
4 th Order gB^T - *j$j$*(B L T) 2 + J^^(^r) 3 


l+§B t T-f(lfeT) a 
S +18**7 ~mm(B L T) 2 


^ + ^ L r + S(^7) a 

i + £BlT - H2&(B L T) 2 £ + f t B L T+^(B L T) 2 


Table 3*8. Loop Filter Constants as Power Series in Bx.T in the Case of Phase and 
Phase-rate Feedback and Computation Delay of One Update Interval 

Super critically damped: rj? = 0, A< = 1, for all roots 

Kl 025 ** Q3S /f? 04 8 If 

l 9t Order IB L T - 24{B L T) 2 + 20&{B L T) S 

Order fB L T - 2»'(£ir) 2 + ‘“'“(BtT) 5 i + i S lT - W^T)* 

B rd Order ff&r- + igig|a ( B I .D» » + gfBtr- {gff}(lW £ + &B t T + igHilBiT)’ g 

4'* Order fBJ- ugm(B L ir + 2§?gJJ{§i?<B,,jr i + gBir-gJfglB,.!-)’ A + 5 fc + ABiT+jgSi(B t r)» 

Standard underdamped: 17? = —1, A» = 1, for all roots 

Ki 

V* Order ABlT - 24{B L T)* + 208(B L T) S 
2 nd Order | Bi r _ W[b l T ) 2 + 2if§2(fli,r) 3 

3rd ° rder ®(^r) a + 5|m2^o(^r) 3 

4** Order |a B lT _ ns2fi (BtT)a + mififp^r) 3 


° >3 H ° ,= ff 

1 + >B t r - ff (B t T)» 

J + gjB^r-JflffflfB.T)’ j, + J§;fli.r + HHg(flir) 1 
$ + ¥,B L T-‘gff(B L T)> i + iB L T-^,(B L Ty X + ± Bl T - 



Table 4-1. Transient-free Initialization of DPLLs 
Based on APriori Information 

For a loop filter specified by update interval T , computation delay n c update intervals, constants K \ , 
K*, Kn and input signal phase known in the form <f>(t) = <£ n +<^(t-t n ) + ...+ 7^0n V) (t-t n ) iV 
at time t n » loop sums and loop model phase 4>m,n are intialized at the end of the “phantom” n th 
interval as indicated. 


V* Order 

= 

2 nd Order 

= 

n-n« 

X) = 
1=1 

3 rd Order 

6 <t> a , = 

n-n„ 

H* = 

t=i 

n-n 0 t 
t=l j=l 




r 3 j(3) 

-^(rvi 2) - km ..) 



Table 4-1. continued Transient-free Initialization of DPLLs 

Based on APriori Information 


4 th Order 

6<f>, a = 

n-rt 0 

= 


T A A<) 

TC A ** 


3^(r 3 ^ 3) - £*L 4) - K z 6<f> aa ) 

J2 J2 6< t>j = ~ “ K 2 6<f> aa - Kzf^ 6<f>i) 

*=1 3 = 1 ,=1 

£ £ £ = i(Wf» + £#<?> + £*<?> + g*S, 4 > - JMfc. - K, £ «* 


t=l 
n-n, > 


*=i y=i fc=i 


«=x 


For all orders: 


{ 4n — Atan extractor 

<t> n - Arcsin(2ir 6<j> aa )/ 2ir, Sin extractor 


*»EE**i 

*=i y=i 
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ABSTRACT OF THE DISCLOSURE 


In a new formulation for digital phase-locked loops, 
loop-filter constants are determined from loop roots 
5 that can each be selectively placed in the s-plane on 
the basis of a new set of parameters, each with simple 
and direct physical meaning in terms of loop noise 
bandwidth, root-specific decay rate, or root-specific 
damping. Loops of first to fourth order are treated in 
10 the continuous-update approximation (B L T + 0) and in a 
discrete-update formulation with arbitrary B L T. 
Deficiencies of the continuous-update approximation in 
large-B L T applications are avoided in the new discrete- 
update formulation. A new method for direct, 

15 transient-free acquisition with third- and fourth-order 
loops can improve the versatility and reliability of 
acquisition with such loops. 
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